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Abstract 

We study the equations obtained from hnearizing the compressible Navier-Stokes 
equations around a steady-state profile with a heavier fluid lying above a lighter fluid 
along a planar interface, i.e. a Rayleigh- Taylor instability. We consider the equations 
with or without surface tension, with the viscosity allowed to depend on the density, 
and in both periodic and non-periodic settings. In the presence of viscosity there 
is no natural variational framework for constructing growing mode solutions to the 
linearized problem. We develop a general method of studying a family of modified 
variational problems in order to produce maximal growing modes. Using these growing 
modes, we construct smooth (when restricted to each fluid domain) solutions to the 
linear equations that grow exponentially in time in Sobolev spaces. We then prove 
an estimate for arbitrary solutions to the linearized equations in terms of the fastest 
possible growth rate for the growing modes. In the periodic setting, we show that 
sufficiently small periodicity avoids instability in the presence of surface tension. 



1 Formulation of the problem 



1.1 Formulation in Eulerian coordinates 

We consider two distinct, immiscible, viscous, compressible, barotropic fluids evolving with 
or without surface tension within the infinite slab f2 := x [—m,i) C M.^ with m,i > 
for time t > 0. The fluids are separated from one another by a moving free boundary 
surface that extends to infinity in every horizontal direction; this surface divides Q into 
two time-dependent, disjoint, open subsets so that Q = ^+{t) U Q_{t) U and 

= n Q-{t). The fluid occupying fl+(t) is called the "upper fluid," and the second 
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fluid, which occupies fl-(t), is called the "lower fluid." The two fluids are described by their 
density and velocity functions, which are given for each t > by 

P±(-, t) : fi±(t) ^ M+ and u±{-, t) : VL±{t) (1.1) 

respectively. We shall assume that at a given time t > the density and velocity functions 
have well-deflned traces onto S(t). 

For t > and x G f2±(t) we require that the fluids satisfy the pair of compressible 
Navier-Stokes equations: 

dtp± + div(p±M±) = ^^2^ 
p±{dtu± + u± ■ Vu±) + div S± = -gp±e3, 

where the viscous stress tensor is given by 

S± = P±{p±)I - £±(p±) (^Du^ + - ^ divM±/^ - <5±(p±)(divM±)J. (1.3) 

In this expression the superscript T means matrix transposition and / is the 3x3 identity 
matrix. The coefficients of viscosity are allowed to vary smoothly with the density, i.e. 
^±y^± £ C°°{{0, oo)), but we assume that the shear viscosity satisfies e± > and that the 
bulk viscosity satisfies S± > 0. In the equations we have written g > for the gravitational 
constant, 63 = (0, 0, 1) for the vertical unit vector, and —ge^ for the gravitational force. We 
have assumed a general pressure law of the form P± = P±{p) > with P± G C°°((0, 00)) 
and strictly increasing. We will also assume that l/P^ ^ Finally, in order to 

create the Rayleigh- Taylor instability, i.e. construct a steady-state solution with an upper 
fiuid of greater density at E(t), we will assume that 

Z ■={ze (0, 00) I P^{z) > P+{z) and P^{z) G P+((0, 00))} ^ 0. (1.4) 

In particular this requires the pressure laws to be distinct, i.e. P_ 7^ P+. For a physical 
discussion of the Rayleigh- Taylor instability, we refer to to [7] and the references therein. 

For two viscous fiuids meeting at a free boundary with surface tension, the standard 
assumptions are that the velocity is continuous across the interface and the jump in the 
normal stress is proportional to the mean curvature of the surface multiplied by the normal 
to the surface (cf. [9]). This requires us to enforce the jump conditions 

(M+)|E(t) - (M-)lE(t) = ^^^^ 

where we have written the normal vector to S(t) as u, and for the trace of a quantity / 
on S(t). Here we take H to be twice the mean curvature of the surface S(t) and the surface 
tension to be a constant a > 0. We will also enforce the no-slip condition at the fixed upper 
and lower boundaries; we implement this via the boundary condition 

M_(xi, X2, —m, t) = u+{xi, X2, i,t) = for all (xi, X2) G M^, t >0. (1.6) 
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The motion of the free interface is coupled to the evolution equations for the fluids (11.21) 
by requiring that the surface be advected with the fluids. More precisely, if V{x, t) G 
denotes the normal velocity of the surface at x G S(t), then V{x, t) = {u{x, t) ■ z/(x, t))z/(x, t), 
where v[x, t) is the unit normal to S(t) at x and u{x, t) is the common trace of u±{-, t) onto 
S(t). These traces agree because of the first jump condition in (11.51) . which also implies that 
there is no possibility of the fluids slipping past each other along S(t). 

To complete the statement of the problem, we must specify initial conditions. We give 
the initial interface S(0) = Eq, which yields the open sets fi±(0) on which we specify the 
initial data for the density and velocity, p±(0) : f2±(0) M"*" and m±(0) : f2±(0) M^, 
respectively. 

It is sometimes desirable to add the additional assumption that solutions are periodic in 
the horizontal directions. More precisely, we can require that for L > 0, the domains fi±(t) 
and the free interface S(t) are horizontally 27rL periodic in that 

f2±(t) = f2±(t) + 2'KLkiei + 27rLA;2e2 and S(t) = S(t) + 2TiLkiei + 2nLk2e2 (1.7) 

for any {ki, k2) G Z^. Then the density and velocity are periodic on Q±{t): 

p±{x + 2'KLkiei + 2-nLk2e2, t) = p±{x, t) for all x G ^±it), (1.8) 

u±{x + 2TiLkiei + 2TrLk2e2, t) = u±{x, t) for all x G Q±{t). (1.9) 



1.2 Reformulation in Lagrangian coordinates 

The movement of the free boundary and the subsequent change of the domains fi±(t) in 
Eulerian coordinates create numerous mathematical difficulties. We circumvent these by 
switching to Lagrangian coordinates so that the interface and the domains stay fixed in time. 
To this end we define the fixed Lagrangian domains Q- = x (— m, 0) and 0+ = x (0, i) 
in the non-periodic case, and Q- = (27rLT)^ x (— m, 0) and fl+ = (27rLT)^ x (0,£) in the 
periodic case. Here we have written 27rLT for the 1— torus of length 27rL. 
We assume that there exist mappings 

: fi± fi±(0) (1.10) 

that are continuous across {x^ = 0}, invertible in the non-periodic case, and invertible on 
their image in the periodic case. We further require that Eq = r]j_{{xs = 0}), '/^^({xa = 
i}) = {xs = i}, and r^^ ({xs = — m}) = {xs = —m}; the first condition means that Eq 
is parameterized by the either of the mappings rj^ restricted to {x^ = 0} (which one is 
irrelevant since they are continuous across the interface), and the latter two conditions mean 
that T]^ map the fixed upper and lower boundaries into themselves. 
Define the fiow maps, t]±, as the solutions to 

dtV±{x, t) = M±(r/±(x, t), t) 
ri{x, 0) = ri^{x). 

We think of the Eulerian coordinates as {y,t) with y = ri{x,t), whereas we think of La- 
grangian coordinates as the fixed {x,t) G f2 x M"*". In order to switch back and forth from 
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Lagrangian to Eulerian coordinates we assume that r)±{-,t) are invertible in the non-periodic 

case and invertible on their images in the periodic case. In the non-periodic case, this imphes 
that ^±(t) = Ti±{Qj-,t), and since and r^^ are all continuous across {x^ = 0}, we have 
S(^) = V±{{^s = 0},t). In other words, the Eulerian domains of upper and lower fluids are 
the image of Q± under the mappings r)± and the free interface is the image of {^3 = 0} under 
the mapping rj±{-,t). In the periodic case, 

n^{t)= y (r^±(0±,t) + 27rLA;iei + 27rLA;2e2), and (1.12) 

(fci,A:2)eZ2 

^{t) = U (^±({^3 = 0}, t) + 2TTLkiei + 2nLk2e2) . (1.13) 

(fci,fc2)ez2 

We define the Lagrangian unknowns 



v±{x,t) = u±{r]±{x,t),t) 
q±{x,t) = p±{v±{x,t),t), 



;i.i4) 



which are defined for {x,t) G fi± x M"*". Since the domains Q± are now fixed, we henceforth 
consolidate notation by writing r],v,q to refer to r}±,v±, q± except when necessary to distin- 
guish the two; when we write an equation for 77, t;, 5 we assume that the equation holds with 
the subscripts added on the domains Define the matrix A via A'^ ~ (Dr])~^, where D 
is the derivative in x coordinates and the superscript T denotes matrix transposition. Then 
in Lagrangian coordinates the evolution equations for v, q, 77 are, writing dj = d/dxj, 

dtTji = Vi 

dtq + qAijdjVi = (1.15) 
qdtVi + AjkdkTij = -gqAijdjTjs, 

where the viscous stress tensor in Lagrangian coordinates, T, is given by 

Tij = P{q)Iij - e{q) (^AjkdkVi + Aikd^Vj - ^(Afe^fe^;/)/^^-^ - 5{q){AikdkVi)Iij. (1.16) 

Here we have written lij for i,j component of the 3x3 identity matrix / and we have 
employed the Einstein convention of summing over repeated indices. 

To write the jump conditions, for a quantity / = /±, we define the interfacial jump as 

I/l :=/+|{.3=o}-/-|{.3=o}- (1-17) 
The jump conditions across the interface are 

W = « (1.18) 



where we have written 

d\ri X d2ri 



n :- 



\di7J X d27]\ 



(1.19) 

{0:3=0} 
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for the unit normal to the surface E(f) = ^{{^a — 0}? ^) H for twice the mean curvature 
of Since S(t) is parameterized by 77, we may employ the standard formula for the mean 
curvature of a parameterized surface to write 



\dir]\ 




2(9i77 • d2'n)did2'n + \d2'nf' dfr] 






2 1 


\92rj\ 


2 


\dir} ■ d2r]\ 
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^= ' " ,2::\2 : ' 2 " -n. (1.20) 



Finally, we require the no-slip boundary condition 

V-{xi,X2, —m, t) — v+{xi,X2, i, t) — 0. (1-21) 

Note that the jump and boundary conditions are the same in the periodic and non-periodic 
cases. 

1.3 Steady-state solution 

We seek a steady-state solution with v = 0,1] = Id, q{x,t) = poix^) with the interface given 
by ri{{x3 = 0}) = {x3 = 0} for all t > 0. Then = 0, n = eg, and A = 7 for all t > 0, and 
the equations reduce to the ODE 

^ - -,Po (1.22) 

dX3 

subject to the jump condition 

[P(po)l = 0. (1.23) 
To solve this we introduce the enthalpy function defined by 

h^{z) = / -^^dr. (1.24) 



The properties of P± guarantee that h± e C°°((0, oo)) are both strictly increasing, and hence 
invertible on their images. The solution to the ODE is then given by 

f ^ _ I h-^i^-iPo) - 9x3), -m<X3<0 
^°^"^-\/.-(Mp^)-^-3), 0<a;3<^. ^ 

where Pq > is a free parameter satisfying P_{pq) e P+((0, 00)), which allows the jump 
condition to be satisfied by choosing > according to 

p+ = p-^(P_(Po))- (1-26) 
For pq to be well-defined on O, we will henceforth assume that £, m > are chosen so that 
(/i_(po ) + gm) e /i_((0, 00)) and (/i+(p+) - gl) e /i+((0, 00)). (1.27) 

Note that po is bounded above and below by positive constants on {—m,i) and that po is 
smooth when restricted to (— m, 0) or (0,£). 
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Since we are interested in Rayleigh- Taylor instability, we want the fluid to be denser 
above the interface, i.e. Pq>Pq. This requires us to choose so that 

^"'(^-(Po )) > Po ^ P-iPo) > PM- (1-28) 

The latter condition is satisfled for any p^ G Z, where Z was deflned by (11.41) : we assume 
Po takes any such value. Then 

Ipo}=P^-Po>0. (1.29) 

For the sake of clarity, we include an example of the solution, pq, when the pressure laws 
correspond to polytropic gas laws, i.e. P±{p) = K^p^^ for K± > 0,7± > 1. The solution is 
then given by 

' [{pov-'-'-^^.y'^'"'^ ^3<o 

\ptV^-'-'-^^^f '''''' < X3 < ^(p+)---^ (1-30) 

.0 ^3>^(Pcr)--^ 

with modiflcation to solutions po(a;3) = p^ ex.\){—gx^/ K±) when either 7+ or 7_ is 1. The 
jump condition requires that 

P^=(^j {PoV-''^- (1-31) 
For a polytropic gas law, the condition that Po > Po is equivalent to 

^) (Po > Po ^ (Po > ^- (1-32) 

If 7+ = 7_ this requires > and any choice of pg > 0. If 7+ 7^ 7- then K^, > 
can be arbitrary, but we must require that pg > satisfles 



Po > XT if 7- > 74 



1/(7+ -7-) 



;i.33) 



Po <{k^) if 7+ > 7-- 

In either case, to avoid the vanishing of po, i is chosen so that 

but the parameter m > may be chosen arbitrarily. 

1.4 Linearization around the steady-state 

We now linearize the equations (ll.lSp around the steady-state solution v = 0, t] = Id, q = po. 
The resulting linearized equations are, writing ri,v,q for the unknowns, 

dtV = V , ^ 

' (1.35) 
dtq + podivi; = 

6 



and 



Podtv + V{P'{po)q) + gqe^ + gpoVr]^ 

= div (^eo (^Dv + Dv^ - ^(divt;)/^ + 6o{diYv)I^ , (1.36) 

where Eq = e{po) and So = S(po). 

The jump conditions hnearize to |f ] = and 

lP'ipo)qI - soiDv + Dv^) - (5o - 2£o/3) div vl} 63 = a A,,,,, 7/363, (1.37) 

while the boundary conditions hnearize to f--(xi, X2, —m, t) = v+{xi, X2,i,t) = 0. We assume 
that initial data are provided as 77(0) = rjo, t>(0) = vq, g(0) = go that satisfy the jump and 
boundary conditions in addition to the assumption that |?7o] = 0, which implies that ri{t) is 
continuous across {0:3 = 0} for all t > 0. 

1.5 Growing mode ansatz 

We will look for a growing normal mode solution to fll.35l) - fll.36l) by first assuming an ansatz 

v{x, t) = w{x)e^\ q{x, t) = q{x)e^\ r]{x, t) = fi{x)e^^ (1.38) 

for some A > 0, which is the same in the upper and lower fluids. Plugging the ansatz into 
fll.35p - fll.36p . we may solve the first and second equations for fj and q in terms of v. Doing 
so and eliminating them from the third equation, we arrive at the time-invariant equation 

X'^Pqw - V(P'(po)po div w) - gpo div ^63 + gpoVw^ 

= div (^Xeo (^Dw + Dw^ - ^(divw)/^ + A5o(divu;)/^ . (1.39) 

This is coupled to the jump conditions {wj = and 

liXSo - 2Xeo/3 + P'(po)po) div wl + XeoiDw + Dw^)] 63 = -aA^.^^^wse^, (1.40) 

and the boundary conditions w_(xi,X2, —m) = w^{xi,X2,i) = 0. Notice that the first jump 
condition implies that the assumptions on ri{0) = fj{0) = w{0)/X mentioned in the last 
section are satisfied. 

Since the coefficients of the linear problem (11.391) only depend on the 0:3 variable, we are 
free to make the further structural assumption that the Xi,X2 dependence of w is given as 
a Fourier mode e*^ where x' ■ ^ = Xi^i + ^2^2 for ^ € in the non-periodic case and 
^ G L~^7j X L~^Z in the 27rL periodic case. Together with the growing mode ansatz, this 
constitutes a "normal mode" ansatz, which is standard in fluid stability analysis [T]. We 
define the new unknowns if,0,ilj : (— m, £) — R according to 

wi{x) = -iip{x3)e'''''^,W2{x) = -i6{x3)e'''''^, and W3{x) = ■?/'(x3)e*''''^. (1-41) 
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The utihty of the new unknowns is seen in the pair of equations 

dww = (ei</^ + 6^ + w^3)e'"'-^ (1.42) 

and 

Dw + Dw^ = I + 2^20 ti^2^-e') \ e"'< (1.43) 

i{i2^-e') 2^' / 

For each fixed ^, and for the new unknowns (/^(xs), 6'(x3), ■{/'(xs), and A we arrive at the 
following system of ODEs (here ' = d/dx^). 

- {\e,ip')' + [AVo + Aeo lei' + il (A5o + A£o/3 + P'(po)po)] V 

= -6 [(A5o + Aeo/3 + P'(po)po) ^' + (A4 - ^7Po)^] - ^16 [A^o + A^o/S + P'(po)po] 

(1.44) 



- (A^o^')' + [AVo + Aeo lei' + ^2 (A5o + A£o/3 + P'(po)po)] ^ 

= -6 [(A5o + Aeo/3 + P'(po)po) ^' + (A4 - 9Po)i^] - [A5o + Aeo/3 + P'(po)po] ^ 

(1.45) 



- [(4Aeo/3 + A5o + P'(po)po) ^']' + (A^o + A^o ^ 

= [(A5o + Aeo/3 + P'(po)po) {^1^ + 6^^)]' + {gpo - KWiV + m (1.46) 

The first jump condition yields jump conditions for the new unknowns: 

M = m = m = 0. (1.47) 

The second jump condition becomes 



(A5o - 2A£o/3 + P'(po)po) {^1^ + ^29 + iP')e, + Xso \ ^(6^ - 9') 



alel'^es, (1.48) 



which implies that 

[A£o(^' - mi = [A^o(^' - mi = (1.49) 

and that 

l{\6o + Aeo/3 + P'(po)po)(^' + ^i^^ + 6^)1 + [A^o (V^' - ^if - 6^)1 = ^ l^l' ^. (1.50) 
The boundary conditions 

ifi-m) = ip{£) = e{-m) = e{e) = tPi-m) = tPii) = (1.51) 

must also hold. 

We can reduce the complexity of the problem by removing the component 6. To do 
this, note that if ip, 9, ifj solve the equations fll.44p - fll.46p for ^ e and A, then for any 
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rotation operator R G S0{2), {(p,9) := R{ip,9) solve the same equations for ^ := i?^ with 
ip, X unchanged. So, by choosing an appropriate rotation, we may assume without loss of 
generality that ^2 = and = |^| > 0. In this setting 9 solves 

-{\eo9J+i\'po + \eo\^\y = 

e{-m) = 9{£) = (1.52) 
91 = [A^o^l = 0. 

Multiplying this equation by 9, integrating over (— m, £), integrating by parts, and using the 
jump conditions then yields 

/ Xeo l^f + (AVo + Aeo Kl')^' = 0, (1.53) 

J —m 

which implies that 9 = since we assume A > 0. This reduces to the pair of equations for 

- X'po^ = -{Xeo^y + lei' (4A£o/3 + X6o + P'(po)po) V 

+ lei [(A5o + A^o/3 + P'(po)po) i'' + (A4 - gpom (1-54) 

- X^poiJ = - [(4A£o/3 + X6o + P'(po)po) + A^o l^l' 

- lei [{{X6o + Xeo/3 + P'(po)po) <^)' + {gpo - A4)v'] (1-55) 

along with the jump conditions 

M = M = [A5o(^'-|e|^)l=0, (1.56) 

[(A5o + Aeo/3 + P'(po)Po)(^' + \i\ ^)1 + [Aeo W - \i\ ^)1 =cy\if ^. (1.57) 
and the boundary conditions 

ip{-m) = ip{i) = ij{-m) = ^{£) = 0. (1.58) 

2 Main results and discussion 

In the absence of viscosity {e = 6 = with modified jump and boundary conditions) and for 
a fixed spatial frequency ^ 7^ 0, the equations fll.54p - fll.55p can be viewed as an eigenvalue 
problem with eigenvalue — A^. Such a problem has a natural variational structure that allows 
for construction of solutions via the direct methods and for a variational characterization of 
the eigenvalue via 

where ^ 

E{^, ^) = l [ P'(po)po(^' + lei - 2gpo lei (2.2) 



2 



m 
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and 



(2.3) 



This variational structure was essential to our analysis in [3], where we showed that A — *■ cxo 
as 1^1 —>■ oo, which led to ill-posedness results for both the inviscid linearized problem and 
the inviscid non-linear problem (equations fll.lSp with e = 6 = 0). 

Unfortunately, when viscosity is present the natural variational structure breaks down 
since A appears quadratically as a multiplier of po and linearly as a multiplier of Eq and 6o 
in fll.54p - fll.55p . This presents no obstacle to a stability analysis once a solution is known 
[Ij since the equations imply a quadratic relationship between A and various integrals of the 
solution, which can be solved for A to determine the sign of 3?A. On the other hand, the 
appearance of A both quadratically and linearly eliminates the capacity to use constrained 
minimization techniques to produce solutions to the equations. 

In order to circumvent this problem and restore the ability to use variational methods, we 
artificially remove the linear dependence on A. To this end, we define the modified viscosities 
e = seq and 6 = s6q, where s > is an arbitrary parameter. We then introduce a family 
{s > 0) of modified problems given by 

- AVo<^ = -(V)' + lel' Ue/3 + 6 + P'(po)po) 



+ 1^1 r6 + E/3 + P'ipo)po)f + is'-gpo)^ (2.4) 



-AW = - 4£73 + 5 + P'(po)poW' +e\^\'^P 



-lel 



6 + Xeo/3 + P'(po)po ) ^ ) + (gpo - e'W 



along with the jump conditions 

M = M = [eV-|ei^)l=0, 

+ e/3 + P'(po)po)(^' + lei ^)] + le (^' - 1^1 
and the boundary conditions 

(p{-m) = <^{i) = tp{-m) = = 



(2.5) 

(2.6) 
(2.7) 

(2.8) 



A solution to the modified problem with X = s corresponds to a solution to the original 
problem. 

Modifying the problem in this way restores the variational structure and allows us to 
apply a constrained minimization to the viscous analog of the energy E defined above (see 
f!3J|) ) to find a solution to fpl4l) - f!23!) with A = A(|^| , s) > when s > is sufficiently small 
and precisely when 



0<|e|<|ele: = 



a 



(2.9) 
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We then further exploit the variational structure to show that A is a continuous function 
and is strictly increasing in s. Using this, we show in Theorem 13.81 that the parameter s can 
be uniquely chosen so that 

s = Xm,s), (2.10) 

which implies that we have found a solution to the original problem fll.54l) -( n..55l) . This 
choice of s allows us to think of A = A(|^|), and gives rise to a solution to the system of 
equations fOD - fOGl) as well. 

Theorem 2.1 (Proved in Section [3A1) . For ^ G so that < |^|^ < g |po] there exists 
a solution cp = (pi^^xs), 6 = e{^,X3), ^ = ^(^^xs), and A = A(|^|) > to (Olj) -( iLiHl ) 
satisfying the appropriate jump and boundary conditions so that ■?/'(^,0) 7^ 0. The solutions 
are smooth when restricted to (— m, 0) or {0,£), and they are equivariant in ^ in the sense 
that if R & 5*0(2) is a rotation operator, then 

/ip{R^,x,)\ /Rn R12 0\ fifitx.y 
0{R^,x,) \ = \R2i R22 e(e,X3) 1 . (2.11) 



U{R^,X3) V 1/ U{^,X3 



Without surface tension (cr = 0) it is possible to construct a solution to fl2.4p - fl2.5p 
with A > for any ^ 7^ 0, but with surface tension (a > 0) there is a critical frequency 
l^lc ~ V 9 [Po] / 0" for which no solution with A > is available if |,^| > |^|^. In the non- 
periodic case, we capture a continuum |^| G (0, |^|^) of growing mode solutions, but in the 
27rL periodic case we only find finitely many. Indeed, if 



-^<L, (2.12) 

then a positive but finite number of spatial frequencies ^ G (L~^Z)^ satisfy |,^| < |^|^. On 
the other hand, if 

then our method fails to construct any growing mode solutions at all. 

It is important to know the behavior of A(|^|) as |.^| varies within < |.^| < |^|^. We show 
in Proposition 13.91 that A(|^|) is continuous and satisfies 

lim A(l^l) = lim Xm) = 0. (2.14) 

In the non-periodic case, this implies that there is a largest growth rate 

0<A:= max A(|^|), (2.15) 

o<l€l<l?L 

and in the periodic case for L satisfying fl2.13p the largest rate is 

0<Al:= sup{A(|e|) I e e (L-'Zf and |^| G (0, (2.16) 

Note that in general < A. In either case, the largest growth rate is achieved for some 
particular choice of ^. 
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The stabihzing effects of viscosity and surface tension are evident in these results. As 
we showed in [1], without viscosity or surface tension, A(|^|) — > oo as |^| oo. With 
viscosity but no surface tension, all spatial frequencies remain unstable, but the growth rate 
A(|(^|) is bounded and decays to as |^| — > oo. With viscosity and surface tension, only 
a critical interval of spatial frequencies are unstable, and A(|^|) remains bounded. Finally, 
with viscosity and surface tension and the periodicity L satisfying (12.131) there do not exist 
any growing modes. 

In the periodic case when L satisfies (12.121) . the solutions to (ll.44p - (ll.46p constructed in 
Theorem 12. II immediately give rise to growing mode solutions to (ll.35l) - (ll.36p . 

Theorem 2.2 (Proved in Section l3.3p . Suppose that L satisfies (12.120 and let ^1,^2 G 
(L^^Z)^ be lattice points such that ^1 = —^2 (I'^d A(|^i|) = Al, where Al is defined by 
(EUD. Define 

wi^, X3) = -iifi^, X3)ei - ie{^, 3:3)62 + 3:3)63, (2.17) 
where ip, 9, ip are the solutions provided by Theorem \2.1\ Writing x' = XiCi + 0:262, we define 

2 

r7(x,^)=6^-*J]^(e„X3)6^^'■«^ (2.18) 
2 

v{x,t) = Ai6^-*5^^&(0,X3)6^^'•^^ (2.19) 
i=i 

and 

2 

q{x,t) = -e^^Vo(a:3)5Z(ei ■e,¥'(e„a;3) + 62 -e.^le^^s) + 53V'(e„^3))e^"'-^^-. (2.20) 

i=i 

Then ri,v,q are real solutions to (11.350 - (11.360 and the corresponding jump and boundary 
conditions. For every t > we have r]{t),v{t),q{t) G H^{Q) and 

\\vmHk = e'''^\\vmH'' (2-21) 

Remark 2.3. In this theorem, the space H''{Q) is not the usual Sobolev space of order k, 
but what we call the piecewise Sobolev space of order k. See (I3.68P for the precise definition. 

In the non-periodic case, although A = A(|^|) for some |^| G (0, |^|^), no L^(fi) solution 
to (ll.35p - (ll.36p may be constructed from a solution to (II. 440 - 01. 461) as in the periodic case 
since 6*^ ^ L'^{Q). We get around this problem by utilizing a Fourier synthesis of such 
solutions. The tradeoff for getting L^(f2) solutions is that the growth rate is not exactly 6^*. 
Nevertheless, it is possible to construct solutions that grow arbitrarily close to this rate. 

Theorem 2.4 (Proved in Section IS75|) . Let f G C^((0, |^|^)) be a real-valued function. For 
^ G with 1^1 G (0, I^IJ define 

u^i^, = -i^i^, 3:3)61 - i9{^, ^3)62 + ipi^, 0:3)63, (2.22) 
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where (f, 6, ip are the solutions provided by Theorem \2.1[ Writing x' = XiCi + X2e2, we define 

r^{x,t) = -^[ /(|e|)^(e,a:3)e^(l«IV^''-«rfe, (2.23) 

v{x,t) = ^ f A(|e|)/(|ei)^(e,X3)e^(l^l)V-'-«^ie, (2.24) 

and 

q(x,t) = -P^ [ fmmM^,^3)+m^,^3) + d.M^,xs))e'^\^\^'e'^'<dC (2.25) 

Then r],v,q are real-valued solutions to the linearized equations (11.351) - fll.36p along with the 
corresponding jump and boundary conditions. The solutions are equivariant in the sense that 
if R & 5*0(3) is a rotation that keeps the vector 63 fixed, then 

rj{Rx, t) = Rri{x, t),v{Rx, t) = Rv{x, t), and q{Rx, t) = q{x, t). (2.26) 

For every k E N we have the estimate 

h(0)|lH^ + ||t;(0)||^. + ||g(0)||^. <cJ [ (1 + 1^1^)'+^ IfiOfd^X^' < 00 (2.27) 



for a constant Ck > depending on the parameters p^, P±, g,a,m,i; moreover, for every 
t > we have r]{t),v{t), q{t) G H'' and 



e*^o(/) h(0)||^. <||r/(t)||^. <e*^h(0)||^. 

e*^«(^Ml^(0)|lH^<||t;(t)L. <e*^||t;(0)||^. (2.28) 



where 



Ao(/) = ,^, inf A(|e|)>0 (2.29) 



and A is given by (12.151) . 



The vertical component of the initial linearized flow map at the interface between the 
two fluids is given in the periodic case by 

r/3(xi, X2, 0, 0) = 2^(6, 0) cos(x' ■ 6), (2.30) 

and in the non-periodic case by 

r/3(xi,X2,0,0) = -^ / fm^{^,0)cos{x'-Odt (2.31) 

Since iIj{^,0) 7^ for any choice of ^, a nonzero / in general gives rise to a nonzero 
773(xi, X2, 0, 0) in the non-periodic case, and ?73(xi, X2, 0, 0) cannot vanish identically in the 
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periodic case. From this we see that vertical displacement is essential to our unstable solu- 
tions. 

It is conceivable that the solutions we construct via the modified viscosity trick somehow 
fail to achieve the fastest growing modes, and so it is not obvious that the growing solutions 
constructed in Theorems 12.21 and 12.41 grow in time at the fastest rate possible. Nevertheless, 
this result is true. In the non-periodic case and in the periodic case when L satisfies f l2.12p . 
we can estimate the growth in time of arbitrary solutions to fll.35p - fll.36p in terms of A and 
A^. The technique we employ was inspired by a similar result, proved in [2J, for the inviscid, 
incompressible regime with smooth density profile. 

To state the result, we first define the weighted norm and the viscosity seminorm by 



Po \v\ 



and 



£o 
n 2 



Dv + Dv'^ - -{divv)I 



+ 5o |divt>| 



(2.32) 



and for i = 1, 2 we write (■, ■)j for the inner-product giving rise to each. 



Theorem 2.5 (Proved in Section . Let v,ri,q be a solution to f ll.35l) -f ri.36p along with 
the corresponding jump and boundary conditions. Then in the non-periodic case 



Ht)\\i + \\v{t)\\i + \\dMml 



< Ce 



2At 



\\dtvml + \\vml + \\vm 



a 



[ \'^.„^MO)f) (2.33) 



for a constant < C = C{pQ,P±,A,e,S,a,g,m,i). In the periodic case with L satisfying 
f l2.12p . the same inequality holds with A replaced with and the integral over replaced 
with an integral over (27rLT)^. 



In the periodic case, when there is surface tension and L satisfies (12. 130 . our method fails 
to construct any growing mode solutions. A priori this does not rule out exponential-in-time 
growth of arbitrary solutions to the linearized equations, but it turns out that exponential 
growth is impossible, and a sort of stability estimate is available. 



Theorem 2.6 (Proved in Section 22]) • In the periodic case let L satisfy f l2.13p . For j > 1 

define the constants Kj > in terms of the initial data via 



Po- 



|9Mo)| 



+ 



^'(Po)po 



div^r t'(o) 
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P'iPo) 

+ 



a 



V.,,.,c^'-^3(0)r. (2.34) 



(27rLT)2 2 



Then solutions to fll.35p - fll.36p satisfy 

imw, + imh < imwi + ii^wii^ + 1 {\\vm, + ii^;(o)iy + 2t'/'^,, (2.35) 

||^^(t)|li + \\vm2 < ll^(0)lli + 11^(0)112 + 3v^v^, (2.36) 
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and for j > 1 

1 ■ ,,2 f°° M • m2 

sup-\\dlv{t)\i+ \\dlv{t)Ldt <2Kj (2.37) 
t>o 2 Jo 

sup||^'t;(t)||' < ||9^(0)||' + 2v^/^. (2.38) 

Our method of studying a family of modified variational problems in order to produce 
growing solutions to linearized problems where viscosity has destroyed the proper variational 
structure is quite general and robust. The method may be used to construct growing so- 
lutions to the compressible Navier- Stokes- Poisson equations with viscosity (cf. [8J), and we 
expect it to be useful for many other viscous, compressible hydrodynamic stability prob- 
lems. The linear instability analysis of this paper comprises the first step in an analysis 
of the non-linear instability of the full equations (11.151) . which will be completed in [5]. A 
non-linear instability analysis of the compressible Navier-Stokes-Poisson equations based on 
linear growing solutions constructed using our method will be completed in [6] for the case 
of constant viscosity and in ^ for the case of density-dependent viscosity. 

The plan of the paper is as follows. In Section 3 we study the family of modified vari- 
ational problems in order to produce growing solutions to (ll.35p - (ll.36p . In Section 4 we 
prove the growth estimates for arbitrary solutions to the linearized problem. 



3 A family of modified variational problems 

3.1 Solutions to (11.441 )— ( [1.45I ) via constrained minimization 



In this section we will produce a solution to (ll.44l) - (ll.45l) with fixed |^| > by first utilizing 
variational methods to construct solutions to the modified problem (I2.4p - (I2.5I) . In order to 
understand A in a variational framework we consider the two energies 



+ \ -m? + {^' -\iwf + +\iwf^ (3.1) 

and ^ 

J{v,i^) = l[ Poiv' + i^'), (3.2) 

which are both well-defined on the space H^{{-mJ)) X H^{{-m,e)). Consider the set 

A={{v,ij) e Hl{{-mA)) X Hl{{-mA)) \ J{v,^) = 1}. (3.3) 

We want to show that the infimum of E{ip, ip) over the set A is achieved and is negative, 
and that the minimizer solves the equations (I2.4l) - (l2.5p along with the corresponding jump 
and boundary conditions. Notice that the jump condition lipj = lipj = holds trivially since 
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(fyip E Hl{{—m^ €)). Also notice that by employing the identity —2ah = (a — hY — (a^ + 6^) 
and the constraint on J(v5, ip) we may rewrite 



E{^,i;) = -g\^\ 



(^(0))^ 



1 

+ 2 



£ ( W - \i\ "^f + W - |{| 4=)' + gl-A' + |{| S')' I > -9 l«l (3.4) 



for any {(p.ip) G ^. Recall that e = se{po), which is smooth when restricted to (— m, 0) 
and (0,£) and bounded above and below by positive quantities for fixed s > 0. In order to 
emphasize the dependence on s G (0, oo) we will sometimes write 



and 



fi{s) := inf E{(p, tp; s) 



(3.5) 
(3.6) 



As the first order of business we show that a minimizer exists. 
Proposition 3.1. E achieves its infimum on A. 

Proof. First note that (13. 4p shows that E is bounded below on A. Let {ipn, ipn) G ^ be a 
minimizing sequence. Then ip^ and tpn are bounded in HQ{{—m,i)) and V'n(O) is bounded 
in M, so up to the extraction of a subsequence (v5n,V'n) ^ ifyi^) weakly in Hq x Hq, and 
{^n,'4'n) — ^ (v^;V^) strougly in x L^. The compact embedding Hq CC H"^^^ ^ C° implies 
that ipn{0) ip{0) as well. Because of the quadratic structure of all the terms in the integrals 
defining E, weak lower semi-continuity and strong convergence imply that 



E{(f, ip) < hm inf E{ipn, "^n) = inf E. 

n^oo A. 

That (v?, ^) G A follows from the strong convergence. 



(3.7) 
□ 



We now show that the minimizer constructed in the previous result satisfies Euler- 
Langrange equations equivalent to fl2.4l) - (l2.5l) . 

Proposition 3.2. Let (</?, ip) & A he the minimizers of E constructed in Proposition \ 3.1[ 
Let fi := E{ip,tlj). Then {^p,ip) are smooth when restricted to (— m, 0) or (0, £) and satisfy 



and 



{e^y + \^fUe/3+~6 + P'{po)po)v 



+ (~6 + e/3 + P'{po)po)^' + {e'-9Po)^ (3.8 



4e/3+~6 + P'{po)po)^' +e\^fi; 



6 + e/3 + P'ipo)po) + (gpo - e')ip 



(3.9) 
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along with the jump conditions 

M = M = [eV-iei^)l = o, 

[{S + e/3 + P'(po)po)(^' + lei V)] + le - iCl ^)1 =ct\^\' V^(0). 
anc? t/ie boundary conditions (f{—m) — ip{£) — ip{—'m) — ip{£) — 0. 
Proof. Fix (¥?o,'0o) e //o((-m,£)) x H^((-m,£)). Define 

r) = J{(p + t(po + r(p,ilj + ti/jo + rip) 
and note that j(0, 0) = 1. Moreover, j is smooth, 



^(0, 0) = y" polv'oV' + V'oV'), and ^(0, 0) = ^ po{ip' 



(3.10) 
(3.11) 

(3.12) 
(3.13) 



So, by the inverse function theorem, wc can solve for r = T{t) in a neighborhood of as a 
function of t so that r(0) = and j{t, r(t)) = 1. We may differentiate the last equation 
to find 

^^-(0,0) + ^(0,0)r'(0) = 0, (3.14) 



dt 



and hence that 



.'(0) = -i|(0.0) = -i 



(3.15) 



Since {(f,ip) are minimizers over A, we may make variations with respect to {(fiQ,ipo) to 
find that 

^ E{ip + tipo + T{t)ip,ilj + tilJo + r{t)ilj), (3.16) 
t=o 



0= -r 
dt 



which implies that 

= a|e|'^(0)(^o(0)+r'(0)^(0)) 

+ f {S + s/3 + P'(po)poW + lel + r'm' + lel 'Po + lel r'(0)^) 

J —m 

g lei Po(V'(<^o + t'(0)(^) + <p{ilJo + t'(O)V')) 
r slV^' - lei (^)(V'; + t'(0)V^' - lel V'o - lei t'(0)(^) 

- lel V^)(¥'o + r'(0)<^' - lel V^o - lel t'(O)V^). (3.17) 



+ 



+ 



Rearranging and plugging in the value of t'(0), we may rewrite this equation as 
a + t {S + e/3 + P'{po)poW + iCl + \C\ Vo) - 9 iCl Poi^Vo + V^o) 



/i /•£ 
e{{^'-\^\<Pm-\^\Vo) + {<p'-\^\i^Wo-\^\^o))^l^ / poi^oV + iPo^) 
■m J —m 

(3.18) 
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where the Lagrange multipher (eigenvalue) is n 
this gives rise to the pair of equations 



= E{ip, ip). Since ipQ and ipQ are independent, 



+ {5 + 4573 + P'(po)po) lei' WO - /' {e lei ^^0)' 

J ~m 



5 + 6/3 + P'(po)Po) ^' + (e^' - fi'Po)^' V^o = / PoWo (3.19) 



and 



cr|e|'V^(0)V^o(0) + 



5 + Ae/3 + P'{po)po )^'+[6 + e/3 + P'(po)po 1^1 V^o 



(e lel V^o)' + / lel' ^ + (^~' - ^?Po) lel v]^o = P Po^^o- (3.20) 



By making variations with y^o, V'o compactly supported in either (— m, 0) or (0, i), we find 
that ip and satisfy the equations fl3.8l) - fl3.9l) in a weak sense in (— m, 0) and (0, £). Standard 
bootstrapping arguments then show that (v^, V^) are in H^{{—m, 0)) (resp. H^{{0,i))) for all 
A; > when restricted to (— m, 0) (resp. (0,£)), and hence the functions are smooth when 
restricted to either interval. This implies that the equations are also classically satisfied on 
(— m, 0) and (0,£). Since G H'^, the traces of the functions and their derivatives are 

well-defined at the endpoints X3 = — m, 0, £. To show that the jump conditions are satisfied 
we make variations with respect to arbitrary (pQ, ipQ G C^{{—m,i)). Integrating the terms 
in (13.191) with derivatives of ipo by parts and using that ip solves (13.81) on (— m, 0) and (0, i), 
we find that 

[£V-|el^)l^o(0) = 0. (3.21) 

Since ipo{0) may be chosen arbitrarily, we deduce the jump condition le{ip' — I^IV')] = 0. 
Performing a similar integration by parts in (I3.20p yields the jump condition 



(3.22) 



The conditions {(pj = {ipj = and ip{—m) = (p{t) = ip{—m) = ip{t) = are satisfied trivially 
since e Hl{{-m,l)) ^ ^^'^^^((-m, £)). 

□ 

We now show that for s sufficiently small, the infimum of E over A is in fact negative. 

Proposition 3.3. Suppose that < |^|^ < g |po] /o"- Then there exists Sq > depending on 
the quantities p^, P±, g,e±, a,m, i, \^\ so that for s < Sq it holds that p{s) < 0. 

Proof. Since both E and J are homogeneous of degree 2 it suffices to show that 

Eiip,^P) 



< 0, 



(3.23) 
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but since J is positive definite, we may reduce to constructing any pair ip) G x Hi 
such that E{(p,ip) < 0. We will assume that (f = |^| so that the first integrand term in 
E{ip,il)) vanishes. We must then construct ip G Hq so that 



(^(0))2 + / gpo^^' + 



(3.24) 

We employ the identity ipip' = {ip'^y/2, an integration by parts, and the fact that po 
solves fll.22p to write 



2l0 



I „/,2 



^?(^(0)) 



gpQ^ 

2 W + f 



Po 



^'(Po) 



(3.25) 



Notice that |po] = Po ~ Po > so that the right hand side is not positive definite. 
For a > 5 we define the test function G HQ{{—m,i)) according to 



Simple calculations then show that 




X3 G [0,£) 

X3 G {—m, 0). 



(3.26) 



v/^(m + £)r(a + l) ^ 
2T{a + 3/2) 

where Oo(l) is a quantity that vanishes as a — >■ 00, and that 

|| + lel^)' + 4(^'^2 



ojl) 



(3.27) 



< C 



(3.28) 



for a constant C depending on a,m,i, \^\. Combining these, we find that 

I A ^ ^ 1^1^ - g [Pol , , 

-E(^a) < h Oq,(1) + sC 



(3.29) 



for a constant C depending on a as well as pQ,P±,g,e±,m,i,\^\. Since a |^|^ < fi'lpo], 
we may then fix a sufficiently large so that the first two terms sum to something strictly 
negative. Then there exists so > depending on the various parameters so that for s < so 
it holds that E{ipa) < 0, thereby proving the result. □ 

Remark 3.4. A simple extension of this argument yields a more quantitative bound that 
holds not only for \^\ fixed, but also uniformly over intervals < a < ^ b < g [po] /a. 
More precisely, there exist two constants Co, Ci > depending on the parameters p^, P±, g, 
e±, a, m, i, a, b so that p(s) < —Cq + sCi for all G [a,b]. 
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The key to the argument presented in this Proposition was constructing a pair ip) so 
that 

^^^^lM(^(0))2 < 0, (3.30) 

which in particular required that ip{0) 7^ and |^|^ < g |po] /cf- We can show that these 
properties are satisfied by the actual minimizers when E{ip,ip) < 0. 

Lemma 3.5. Suppose that {ip,ip) G A satisfy E{(p,tlj) < 0. Then ip{0) 7^ and < 
9 [Pol /f^- 

Proof. A completion of the square allows us to write 
P'{po)po{f+\^\^Y-2gpo\^\^^ 

V^nPo)io{^' + \^\v)-^^,^] +2^/Po^^'--|^^^. (3.31) 
Integrating by parts as in (13.251) . we know that 

29Poi^^' - Sf^^P' = -9 [Pol (^(0))^ (3.32) 
Combining these equalities, we can rewrite E{ip, tp) as 
1 



+ 2 j Ji.i.v'-mY + {i^'-\i\vY)+ 2 WO)) - (3-33) 

From the non-negativity of the integrals, we deduce that if -E(v?, ^) < 0, then ^/'(O) 7^ and 

\i?<9lPoll^. □ 

The next result establishes continuity and monotonicity properties of the eigenvalue p{s). 
Proposition 3.6. Let p : (0, 00) ^ M 6e given by i \3.6\} . Then the following hold. 

1. p E C^^{{^, 00)); o.nd in particular p G C°((0, 00)). 

2. There exists a positive constant C2 = C2{p^, P±, g,£±, o',m, £) so that 

p{s)>-g\^\ + sC2. (3.34) 

3. p{s) is strictly increasing. 

Proof. Fix a compact interval Q = [a,b] CC (0, 00), and fix any pair [ipo, ipo) G A. We may 
decompose E according to 

E{ip,^lj;s)=Eo{f,^P)+sEi{^,i;) (3.35) 
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for 

a 

Kr^l in. ih) := — , ^ , 

2 " 2 

and 



Eo(¥^,^):=^(^(0))' + ^ / P'(po)po(V^'+Kk)'-2^|e|po¥^^ (3.36) 



E,{^,^):=\ [ {5o + eo/3W+\^\^r + eo{{^'-\^\^Pr + i^'-\^\^Y)>0. (3.37) 

^ J -m 

The non-negativity of -Ei implies that E is non-decreasing in s with (yj, ?/)) G .4. kept fixed. 
Now, by Proposition 13. ![ for each s G (0, oo) we can find a pair {{ps,ips) G ^ so that 

E{ips,i^s;s) = inf E{ip,ij;s)=fi{s). (3.38) 

We deduce from the non-negativity of Ei, the minimahty of {ipsyips), and the equahty (13.40 
that 

E(v?o, ^o; ^) > ^(v^o, ^o; s) > ^(v^s, ips] s) > sEi{ips, ips) - 9 1^1 (3.39) 

for all s E Q. This implies that there exists a constant Q < K = K{a,b,(fo,'ipo, 9, |^|) < oo 
so that 

snp E^{ips,^IJs)<K. (3.40) 

Let Sj G Q for i = 1,2. Using the minimality of {ipsi,ipsi) compared to (v^^j, ^/^^J, we 
know that 

/i(si) = E{ips,,'^si;si) < E{ips2,ilJs,;si), (3.41) 
but from our decomposition fl3.35p . we may bound 

E{'fs2,'i/js2;Sl) < E{ips2,1ps2]S2) + |Si - S2\Ei{ips2,^s2) 

= Ks2) + \si - S2I Ei{Lps2,ips2)- (3.42) 
Chaining these two inequalities together and employing fl3.40p . we find that 

f^i-si) < fi{s2) + K \si - S2\ . (3.43) 

Reversing the role of the indices 1 and 2 in the derivation of this inequality gives the same 
bound with the indices switched. We deduce that 

|/i(si)-/x(s2)| <K|si-S2|, (3.44) 

which proves the first assertion. 

To prove (13.340 we note that equality (13.40 and the non-negativity of Ei imply that 

Ms) >-9\^\ + s^ inf E,{<f,^). (3.45) 

It is a simple matter to see that this infimum, which we call the constant C2, is positive. 
Finally, to prove the third assertion, note that if < si < S2 < cxd, then the decomposition 
(I3.35P implies that 

H{si) = E(v5si,^^i; si) < E{(ps2,^s2; si) < E{lps2,^ps2■, S2) = Ai(s2). (3.46) 
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This shows that fi is non-decreasing in s. Now suppose by way of contradiction that = 
/i(s2). Then the previous inequality implies that 

SiEi{ips2,'^s2) = S2Ei{ips2,'ips2), (3.47) 

which means that Ei{(fs2, "^52) = 0- This in turn forces ips2 = V's2 = O5 which contradicts the 
fact that {ips2,ips2) ^ Hence equahty cannot be achieved, and /i is strictly increasing in 
s. 

□ 

Now we know that when < |^|^ < g |po] /o", the eigenvalue /x(s) is a continuous function. 
We can then define the open set 

5 = /i-^((-oo,0)) C (0,00), (3.48) 

on which we can calculate A = ^— /i > 0. Note that S is non-empty by Proposition 13.31 

We can now state a result giving the existence of solutions to f l2.4p -( l2r5l) for these values 
of 1^1 , s. To emphasize the dependence on the parameters, we write 

^ = ^sm,X3),^ = tpsm,X3), and A = A(|e|,s). (3.49) 

Proposition 3.7. For each s G 5 and < |^|^ < g |pol /c^ there exists a solution ^Ps{\^\ , X3), 
"^sd^l ^x^) with A = A(|^| ,s)>Oto the problem fl2.4p - fl2.5l) along with the corresponding 
jump and boundary conditions. For these solutions 'ips{\C,\,0) 7^ and the solutions are 
smooth when restricted to either (— m,0) or (0,£). 

Proof. Let (v^sd^l ,-),4's{\^\ ,■)) ^ be the solutions to f l3.8p -( l3^ constructed in Propo- 
sition 13. 2[ Since s G iS we may write p = — A^ for A > 0, which means that the pair 
((/?,(|^| , ■), V'.dCl , ■)) solve the problem (^-(^. The fact that Vs(KI , 0) 7^ follows from 
Lemma 13.51 □ 

In order for these solutions to give rise to solutions to the original problem, we must be 
able to find s G 5 so that s = A(|^| ,s). It turns out that the set S is sufficiently large to 
accomplish this. 

Theorem 3.8. There exists a unique s & S so that X{\^\ , s) = a/— plsj > and 

s = Xm,s). (3.50) 

Proof. According to Remark l3.4l we know that fi(s] < —Cq+sCi. Moreover, the lower bound 
fl3.34p implies that /i(s) +00 as s ^ 00. Since p is continuous and strictly increasing, 
there exists s* G (0, 00) so that 

5 = /i-^((-oo,0)) = (0,s,). (3.51) 

Since p < on 5, we may define A = y/—^ there. Now define the function $ : (0, s*) 
(0, 00) according to 

^s)=s/m\,s). (3.52) 

It is a simple matter to check that the continuity and monotonicity properties of /i are 
inherited by $, i.e. $ is continuous and strictly increasing in s. Also, lim^^o ^(■s) = and 
lim^^s, $(s) = +00. Then by the intermediate value theorem, there exists s G (0, s*) so that 
$(s) = 1, i.e. s = A(|^| , s). This s is unique since $ is strictly increasing. □ 
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We may now use Theorem l3.8l to think of s = since for each fixed < < g |pol /c 
we can uniquely find s G 5 so that fl3.50p holds. As such we may also write A = A(|^|) from 
now on. 

Using this new notation and the solutions to fl2.4p - fl2.5l) given by Proposition 13.71 we can 
construct solutions to the system system fll.44l) - fll.46l) as well. 

Proof of Theorem \2.1[ We may find a rotation operator R G 5*0(2) so that RC, = (|^| ,0). 
For s = s(l^l) given by Theorem 13.81 define ((/)(^, X3), ^(^, xs)) = , X3), 0) and 

i^i^y^s) = V'sd^l 5^3); where the functions ^9^(1^1 5X3) and i'sil^l ,^3) are the solutions from 
Proposition 13.71 This gives a solution to fll.44p - fll.46p . The equivariance in ^ follows from 
the definition. □ 



3.2 Behavior of the solutions with respect to ^ 

In this section we shall study the behavior of the solutions from Theorem 12. II in terms of ^. 
We assume throughout that |^| G (0, with |^|^ = a/^IpoJT^- The results are primarily 
needed in the non-periodic case, when there is a continuum of spatial frequencies in (0, |^|^). 
The first result shows that A is a bounded, continuous function of |^|. 

Proposition 3.9. The function A : (0, \C,\^) —>■ {0,oo) is bounded, continuous, and satisfies 

lim A(lel) = lim Xm) =0. (3.53) 

Proof. We begin by proving the continuity claim. Since A = y/^JI it suffices to prove the 
continuity of fi = /i(|^|). By Proposition 13.21 for every |^| G (0, there exist functions 
(V'l?!,^!?!) e A satisfying ([3l])-([3lD so that /i(|^|) = E((/?|5|, ^i^i). We have that /i(|^|) < 0, 
which, when combined with (13.40 . yields the bound 

-9\^\ + f f ((^1 - lel %|)' + (V^isi - mf) < /^(lel) < (3.54) 

for all |^|. 

Now suppose G (0, g |pol /o") is a sequence so that |^|^ — >• |^| G (0, g |pol /o")- We may 
assume without loss of generality that I^L e [1^1 /2, (|el + ieL)/2] if a > Oor |e|„ G /2,2 
if cr = 0. In order to make use of the bound (13.540 we must show that s(|^|„) is bounded 
uniformly from below as n ^ 00. By Remark 13.41 there exist positive constants Co, C\ so 
that pdeij < -Co + .(lelJC^i, but -pdeij = A^delJ = ^'(leij, so 

0<.^(|e|J+Ci.(|e|J-Co (3.55) 

and hence s(|^|^) is bounded below by a positive constant. Then (13.540 and the fact that 
((y9|^l ) G A imply that and '^\^\ are uniformly bounded in H^{{—m,£)). Plugging 

into the ODE fl3.8p - fl3.9p in the intervals (— m, 0) and (0, i) separately, we find that and 
w\^\ are uniformly bounded in iJ^((— m,0)) and if^((0,£)). So, up to the extraction of a 
subsequence we have that 

(¥.151^,^^151 J ^ im^^m) strongly in H\i-m,0)) and H\iOJ)). (3.56) 
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This implies that along the subsequence 

MlelJ = ^1^1 J - ^i^i) = Miel). (3.57) 

Since this must hold for any such extracted subsequence, we deduce that A^d^l) 
for the original sequence |^|^ as well, and hence fi is continuous. 

We now derive the limits as |^| 0, |^|^. By fl3.54p . < A^(|,^|) < g\^\ , which establishes 
that limi^i^o A(l^l) = 0. By ([32]) we know that 

(%|(0))' < (3.58) 

but by fl3.33p we also know that 



|2 



A^(|e|)<^M_^(^,^,(0))^ (3.59) 



Chaining the two inequalities together then shows that lim|g|^|^| A(|^|) = 0. 



□ 



Remark 3.10. A trivial consequence of this result is that the supremum of X is achieved. 
We denote the supremum in the non-periodic case by A (see (12.151) ) and in the periodic case 
byAL (see (^M)- 

The next result provides an estimate for the H'^ norm of the solutions {(p, 6, ip) constructed 
in Theorem 12. H which will be useful later when such solutions are integrated in a Fourier 
synthesis. 

Lemma 3.11. Suppose < a < b < |^|^ and that |^| G [a,b]. Let {(f,6,ilj) be the solutions 
constructed in Theorem \2.1\ Then for each k > there exists a constant Ak > depending 
on the parameters a,b, p"^, P±, g,e±,6±, a, m, i, 

M^,-)\\ + \m.-)\\ + \m.-)\\ 

+ Mi.-)\\HHm) + ll^(^'-)ll//^((o/)) + II V'(e,-) 11/^^(0/)) < ^k- (3.60) 

Also, there exists a Bq > depending on the same parameters so that 

Vv'i^r) + 0\^,-)+n^,-) ^^^^ > Bo. (3.61) 

Proof. Since the solutions in Theorem 12.11 are constructed from rotations of the solutions 
constructed in Proposition 13.71 it suffices to prove 



ll^del 

+ Ikdei , -nHHm) + 11^(1^1 ' -nHHm) ^ (3-62) 

for the solutions (p = (p{\^\ ,2:3), ip = p>{\^\ ,0:3) constructed in the theorem. For simphcity 
we will prove an estimate of the H'' norms only on the interval (0, i). A bound on (— m, 0) 
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follows similarly, and the result follows by adding the two. Recall that po cind -P'(po) are 
smooth on each interval (0,£) and (— m,0) and bounded above and below. 

We proceed by induction on k. For k = the fact that (v^d^l , > ")) ^ ^ implies 

that there is a constant Aq > depending on the various parameters so that 

11^(1^1 , -nmrn) + , ■)\\LHm) < ^o- (3.63) 

Suppose now that the bound holds some k > 0, i.e. 

, -nHHm) + ll^d^l ' -^Wnnm) ^ ^'^^ (3-64) 

By Proposition 13. 9[ A(|^|) = s(|,^|) is bounded above and below by positive quantities as 
functions of |^|. Then by differentiating the equations fl2.4p - fl2.5l) we have that there exists 
a constant C > depending on the various parameters so that 

ll^(lei,-)|| + 11^(1^1, Oil 

< C{ym , OIIhhco/)) + > OIlH^ao/))) < CAk := Ak+,. (3.65) 

Then the bound holds for k + 1, and so by induction the bound holds for all k > 0. 

To prove (13.611) we again utilize the fact that \ , ■), ■j/'d^l , ■)) G A. Since po is bounded 
above and below, the bound follows. 

□ 



3.3 Solutions to (|L35D-(fL36D 

In this section we will construct growing solutions to (ll.35p - (ll.36l) by using the solutions to 
(ll.44p - fll.46p constructed in Theorem 12. 1[ In the periodic case this can only be done when L 
satisfies (12.120 . but the construction is essentially trivial since normal mode solutions are in 
L^(fi). In the non-periodic case, we must resort to a Fourier synthesis of the normal modes 
in order to produce L'^{Q) solutions. 

We begin by defining some terms. For a function / G L^(f2), we define the horizontal 
Fourier transform in the non-periodic case via 

f{^i,^2,x,)= [ fix,,X2,Xs)e-'^^'^^+^'^''^dx,dx2 (3.66) 

for ,^ G M^. In the periodic case the integral over must be replaced with an integral over 
(27rLT)^ for ^ G (L~^Z)^. In the non-periodic case, by the Fubini and Parseval theorems, 
we have that / G L^(fi) and 

f \f{x)fdx=^ I |/(e,X3)|'rferfx3. (3.67) 

The periodic case replaces 47r^ with 47r^L^ and the integral with a sum over (L^^Z)^ on the 
right hand side. 
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We now define the piecewise Sobolev spaces. For a function / defined on f2 we write /+ 
for tlie restriction to and /_ for the restriction to For A; G N, define the piecewise 
Sobolev space of order k by 

H\n) = {f\f+e H'in+), /_ G H\n^)} (3.68) 

endowed with the norm = + \\f\\H-{n^)- Writing /_ = (-m, 0) and /+ = 

(0, i), we can take the norms to be given as 



k 

E / (l + l^l')'"'' diM^r) ' (3.69) 



k 

H''{n±) •- 

j=0 



in the non-periodic case; for the periodic case we replace the integral over with a sum 
over (L~^Z)^ on the right hand side. The main difference between the piecewise Sobolev 
space H^{Q) and the usual Sobolev space is that we do not require functions in the piecewise 
space to have weak derivatives across the interface {x^ = 0}. 

The 2itL periodic growing mode solutions may now be constructed. 



Proof of Theorem \2.^ It is clear that rj, v, q defined in this way are solutions to fll.35p -(ITT 



That they are real- valued follows from the equivariance in ^ stated in Theorem 12. 1[ The 
solutions are in H'^IQ) at t = because of Lemma [3. Ill The growth in time stated in (12.211) 
follows from the definition of rj, v, q. □ 

In the non-periodic case the exponentials e*^ '^ are not in L^(n), so we must utilize a 
Fourier synthesis. The tradeoff for utilizing such a synthesis is that the growth rate is not 
exactly e^*, but can be made arbitrarily close to it. 



Proof of Theorem\2^ For each fixed ^ G so that |^| G (0, 

r/(a;,t) = /(|e|)t&(e,a;3)e'(l«l)*e«'-«, (3.70) 

v{x, t) = A(|e|)/(|el)^(e, X3)e^(l«IV-''-«, and (3.71) 
q{x,t) = -po(x3)/(|el)(6<^(e, 2:3) +6^(^,2:3) + 53^(e,a;3))e'(l«l)V^'-« (3.72) 
constitute a solution to f ll.35p - p.36p . Since supp(/) CC (0, |^|^), Lemma 13 . 1 II implies that 

sup \\d^^w{^, ■)\\^^ < 00 for all keN. (3.73) 
Cesupp(/) 

These bounds, the definition of A, and the dominated convergence theorem imply that the 
Fourier synthesis of these solutions given by fl2.23p -( l2.25p is also a solution that is smooth 
when restricted to Q±. The Fourier synthesis is real- valued because /(|^|) is real- valued and 
radial and because of the equivariance in ^ given in Theorem 12. 1[ This equivariance in ^ also 
implies the equivariance of r],v,q written in (I2.26p . 

The bound (I2.27P follows by applying Lemma [3. Ill with arbitrary A; > and utilizing the 
fact that / is compactly supported. The compact support of / also implies that Ao(/) > 0, 
so that Ao(/) < A(l^l) < A for |^| G supp(/). This then yields the bounds (I22H])- □ 
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4 Growth of solutions to the linearized problem 
4.1 Preliminary estimates 

In this section we will prove estimates for the growth in time of arbitrary solutions to 
fll.35p - fll.36p in terms of the largest growing mode: A in the non-periodic case and A.^ in the 
periodic case, defined by (12. 150 and (12.160 respectively. To this end, we suppose that 77, f , q 
are real-valued solutions to f 11.35p - fll.36p along with the corresponding jump and boundary 
conditions (of course, by linearity, we may also handle complex solutions by taking the real 
and complex parts and proceeding with an analysis of each part). 

It will be convenient to work with a second-order formulation of the equations. To arrive 
at this, we differentiate the third equation in time and eliminate the q and r] terms using the 
other equations. This yields the equation 



p^dttv - V{P\po)po divv) + fi'PoVfs - gpodivves 



div ( eo { Ddtv + Ddtv^ - -{div dtv)I ] + 5oidwdtv)I ] (4.1) 



coupled to the jump conditions 
and 



Idtvj = 



(4.2) 



P'(po)po divv)I + EoiDdtV + Ddtv'^) + (5o - 2£o/3) divS^t;/] 63 = -a^,,^,,v^e^. (4.3) 



The function dtv also satisfies dtv{xi^X2^ —m, t) = dtv{xi, X2, i,t) = at the upper and lower 
boundaries. The initial data for dtv{0) is given in terms of the initial data g(0),f(0), and 
77(0) via the third linear equation, i.e. dtv{0) satisfies 

podtv{0) = -gqiO)es - gpo'^VsiO) 



+ div i^eo i^Dv{0) + Dv{Of - -(divf (0))/j + 5o{divv{0))Ij . (4.4) 

Our first result gives an energy and its evolution equation for solutions to the second-order 
problem. 

Lemma 4.1. Let v solve (14. ip and the corresponding jump and boundary conditions. Then 
in the non-periodic case, 



d, 



Po 



\dtvf , P'(po)po 



+ 



div V — 



P'iPo) 



V3 



+ 



+ / 6o\divdtv\ 



£0 



Ddtv + Ddtv^ - -{div dtv)I 



1^3 1 



-^\^xux2n\ ■ (4.5) 



In the periodic case, the same equation holds with the integral over replaced with an 
integral over (27rLT)^. 
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Proof. We will prove the result in the non-periodic case. The periodic case follows similarly. 
Recall that fl+ = x (0,£). Take the dot product of fl4.ip with dtv{t) and integrate over 
After integrating by parts and utilizing (11.221) . we get 

r 2 

/ podtv ■ dttv + P'(po)po(div v){div dtv) - gpo{v3 div dtv + dm div v) + --j^v^dtv^ 
Jn+ P{po) 



+ 



£o 
n+ 2 



+ / 5o|div9tf| 



dtv^- f P[ip+)p+diYvdtVs- [ Tes-dtv (4.6) 

JM2 JK2 



where we have written 



T = (P'(po)po div v)I + eo (^Ddtv + Ddtv'^ - ^ (div dtv)?j + 5^ div dtvl. (4.7) 

We may pull time derivatives out of the first integrals on each side of the equation to arrive 
at the equality 



/ po— h 



\dtv\^ , P'(po)po 



div V — 



9 



+ 1 ^ 

n+ 2 



P'iPo) 

+ 6o \divdtvf = dt gp^ 
Jn+ JR2 



Ddtv + Ddtv'^ - -{diYdtv)! 



I ^3 I 



- / Tes ■ dtV. {Al 



A similar result holds on f2_ = x (— m, 0) with the opposite sign on the right hand side. 
Adding the two together yields 



o , ^ \dtv\ P'ipo)po 



div V 



9 



P'iPo) 



V3 



+ I ^ 

n 2 



Ddtv + Ddtv'^ - -{divdtv)! 



+ I 6o \diYdtvf = dt g |pol 



lTe;-dtvj. (4.9) 



Using the jump conditions, we find that 



|Te3 ■ dtvj = / crA^^^^^vsdtVs 
Jm.2 

= -0- / Vx^,x2V2,-Vx^,x2dtVz = -dt I 7:\^xi,x2Vz? ■ (4-10) 
Jr2 Jir2 Z 

The result follows by plugging this in above. □ 

The next result allows us to estimate the energy in terms of A, which was given by f l2.15p . 
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Lemma 4.2. Let v G H^{Q) be so that v{xi,X2, —m) = v{xi,X2,i) = 0. In the non-periodic 
case we have the inequality 



9 [Pol I |2 



2 ' ^' 2 

< 



IV7 \^ 
A2 



P'iPo)po 



n 



divf 



P'iPo 



-V3 



+2^2 



Dv + Dv'^ - -{divv)I 



+ 6o\dwvf. (4.11) 



In the periodic case, if ^ cy ji^g |po]) < ^, then the same inequality holds with the integral 
replaced with an integral over (27rLT)^ and A replaced with A^. 

Proof. We will again prove only the non-periodic version. Take the horizontal Fourier trans- 
form and apply fl3.67p to see that 



9 Ipo 



V3 



P'iPo)Po 



divi; 



9 



P'iPo) 



V3 



9lPo}-^\^\\ . ,2 
F3 



7, rsl 



P'iPo)po 



P{po) 



t TDI 



P'iPo)po 



i^m + i^2V2 + dsvs 



P'iPu 



-f3 



d^dxs 

dx^ I d^. 
(4.12) 



Consider now the last integrand for fixed ^ 7^ 0, writing (p^x^) = ivi{C,, X3), O^x^) = iv2{C,, X3) 
ipi^xs) = V3{^,X3). That is, define 



9lPo}-cT\^\\ ,,2 



1^1- 



e TDI 



P'{Po)po 



^1^ + ^20 + 4^'--^^ 

P{po) 



where ' = 83. By sphtting 



dx3 (4.13) 



(4.14) 



we may reduce to bounding Z when yj, 6', il) are real- valued functions, and then apply the 
bound to the real and imaginary parts of 6', 

The expression for Z is invariant under simultaneous rotations of ^ and ((y9, Q\ so without 
loss of generality we may assume that ^ = (|^| , 0) with |^| > and 6^ = 0. If a > then we 
assume for now that |,^| < |^|^ as well. Then, using (13. ip with e = A(|,^|)£o and 6 = A(|,^|)5o, 
we may rewrite 



Z{^,9,^p■,0 = -E{^,^■,Km■ + 



So\^'+\^\^\ 



+ 



1 



eo ( iv^' - lei + - lekr + 3 1^' + lei ) (4.15) 
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and hence 



^ J —m ^ J —m \ 



mv + w-\^wv + -^w + \iwv] (4.16) 



For \^\ > the expression for Z is non-positive, so the previous inequahty holds trivially, 
and so we deduce that it holds for all |^| > 0. 

Translating the inequality back to the original notation for fixed ^, we find 



9lP0l-CT\C\\ . ,2 
7^ "^3 



P'{Po)po 



A 



A 



P'iPo) 



V3 



< TT / Pol^^l +77/ So\'i'^lVl+i^2V2 + d3V3\ + — 



|2 _|_ ^0 


B 


' 2 





where 



B = Dv + Dv^ - -{divv)I. 



(4.17) 



(4.18) 



Integrating each side of this inequality over all ^ G and using fl3.67p then proves the 
result. □ 

When 0" > and L is sufficiently small, a better result is available in the periodic case. 

Lemma 4.3. Let v G H^{Q) be so that v{xi,X2, —m) = v{xi,X2,i) = and suppose in the 
periodic case that L satisfies (12.131) . Then 



9 Ipo 



IV7 \^ 



(27rLT)2 



^'(Po)po 



divf — 



P'{Po 



< 0. 



(4.19) 



Proof. Apply the horizontal Fourier transform to see 



ajpo 

(27rLT)2 2 



a 



^'(Po)po 



divf — 



P'{Po) 



ee(L-iz)2 



9lp4-<^\i\\ ^ ,2 

7. 



^'(Po)po 



iiivi + ii2V2 + d^vj, 



P'iPo) 



V3 



dx3. (4.20) 



Because of (127[3|) . the only ^ G {L-^Z^ for which g |pol - 9 > is ^ = 0. Since all but 
the ^ = term on the right side of the last equation are non-positive, we reduce to showing 
that 



9 [Pol |. |2 



^'(Po)po 



d3V3 



9 



P'iPo 



-V3 



dx3 < 0. 



(4.21) 
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For this we expand the term in the integral and integrate by parts to get 



glpol |. |2 



9 



dx-\ = 



which yields the desired inequality. 



P'ipo)po\d3Vs\\ (4.22) 



□ 



4.2 Proof of Theorems [23] and [2761 

With the preliminary estimates in place, we can now prove bounds for the growth of arbitrary 
solutions to (14.1 p in terms of A and A^. 



Proof of Theorem \2.5[ Again, we will only prove the non-periodic case. Integrate the result 
of Lemma 14.11 in time from to t to find that 



\dMt)\' , 
Po^— + 



<Kn + 



Jn 

glPol 



£o 
2 



Ddtv{s) + Ddtv{sf - -{div dtv{s))I 



+ So \div dtv{s)\ ds 



2 Mt)\^ -^\^x,,X2V3{t)\^ 



P'{Po)po 



divt>(t) 



P'{p^ 



Mt) 



(4.23) 



where 
Ko = 



p ftvm ^ I p'ipo)po 



divt;(0) 



P'iPo 



+ I f |V.„.,t;3(0)|^ (4.24) 



We may then apply Lemma 14.21 to get the inequality 

\dMt)f 



Po- ^ 
n ^ 



/ / ? {^^'""^'^ + Ddtv{sf - '^{diY dtv{s))I^ + 6o \diY dtv{s)fds 



+ 



A f So 



2 In 2 



<Ko + ^ I PoHt)]'' 
n 



Ddtv{t) + Ddtvity - -{div dtv it)) I + 5q |divt;(t)r . (4.25) 



Using the definitions of the norms 
previous inequality as 



1 



i,\\dtvmi 



^ , II -112 given in fl2.32p . we may compactly rewrite the 
A2 „ A 

(4.26) 



2 • 



Integrating in time and using Cauchy's inequality, we may bound 



A||i;(t)||^=A||^;(0)||^ + A / 2{v{s),dtv{s))2ds 



t ft 

2 



<A||t;(0)||^+ / \\dtv{s)\\lds + A' \\v{s%ds. (4.27) 
Jo Jo 
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On the other hand 

Kdt \\v{t)\\\ = A2{dtv{t),v{t)), < A' \\v{t)\\l + WdMml ■ (4.28) 
We may combine these two inequahties with (14.261) to derive the differential inequahty 

dt\\v{t)\\l + \\v{t)\\l<Ki + 2A\\v{t)\\l + 2A[ \\v{s)\\lds (4.29) 

Jo 

for Ki = 2Kq/ A + 2 ||f (0)||2 . An apphcation of Gronwall then shows that 



(4.30) 



for all t > 0. To derive the corresponding bound for ||f(t)||2 and ||c?ti'(t)||^ we return to 
(1426]) and plug in <KTj\ and fOOjl to see that 

\ \\dtvm\ + \Ht)\\l < K, + A ||^;(t)||; + 2A /* Hs)\\lds < e'^' (2A ||^;(0)||? + K,) . (4.31) 
The result follows by noting that 







Ko<c [\\dMom + Hmt + Hmt + ^ J^^ i'^.u^mwj (4.32) 

for a constant C > depending on p^, P±, A, e±, 6±, a, g,m, i. 

□ 

In the periodic case when L satisfies (12.13^ we may use Lemma IT3] to improve the above 
result. 

Proof of Theorem \2.6{ . We again integrate the result of Lemma 14.11 in time from to t to 
find that 



Ddtv{s) + Ddtv{sf - -{dwdtv{s))I 

3 



n ^ Jo Jn 

^ , f ^[PoI| .,x|2 1^ .,x|2 f P'{P0)P0 

J(27rLT)2 ^ ^ Jo, 



2 



+ 5q |div9tf (s)I^ ds 



2 



2 



(4.33) 



We may apply Lemma 14.31 to see that all of the integrals on the right side of the previous 
inequality are non-positive, and hence 

\\\dtv{t)\\\ + j\\d,v{s)\\lds<K,, (4.34) 

where the norms are defined by (12.321) . From this we deduce that 

\\v{t)\\, + \\v{t)\\, < \\vmi + 11^(0)112 + 3v^v^. (4.35) 
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Then, using that dtt] = v, we get 

ii + ii^wii2 < ii^(o)iii + imh+tiWvm, + ii^(o)iy + 2t^/'v^. (4.36) 



To derive the estimates for dfV for j > 2 we apply to (14.11) . Then w = d^v satisfies the 
same equation and boundary conditions as v, which allows us to argue as above to derive 
the inequality 

l\\dMt)\\l + 1^ \\dMs)\\lds < (4.37) 
for all j > I. This trivially implies (12.371) . To get (12.381) we bound 



|2 , II ii2 



\div{t)\\l<\\d'tv{0)\\; + 2 I \\div{s)\\Jdi+'v{s)\\^ds 







<\\dM0)\\l + 2(^l^\\dlv{s)\\lds^ ' (^1^ \\dt'v{s)\\lds^ ' 



< ||a^i;(0) 11^ + 2 v/i?;/??^. (4.38) 

□ 
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